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Abstract: We investigate the complete next-to-leading order (NLO) QCD corrections 
to T(1S) + W production at the LHC, and predict theoretically the distribution of the 
T(1S) transverse momentum. We analyse the contributions from different components up 
to the QCD NLO in the pp — > T(1S) +W + X process. Our results show that the pp — > 
T(1S) + W + X process has a large production rate and could be potentially detected at the 
LHC. We see that the differential cross section for the T(1S*) direct-production at the LO 
is significantly enhanced by the QCD corrections, and the 66[ 3 5^] contribution component 

Yds) 

dominates in the whole plotted p T region. We have also calculated the T(1S) meson 
indirect-productions via feed-down decays of T(2S), T(35), \biiXP), Xb2(lP), Xbi(2P), 
and Xb2(2P) mesons. We find that the T(IS') indirect-productions can give important 
contributions to the distribution of p T for the pp — > T(15) + W + X process at the 
NLO. We conclude that the studying the T(1S) + W production at the LHC could provide 
an interesting opportunity in testing the nonrelativistic QCD factorization formalism. 
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1. Introduction 

The study of heavy quarkonium is one of the interesting subjects in both theoretical and 
experimental physics, which offers a good testing ground for investigating the Quantum 
Chromodynamics (QCD) in both the perturbative and non-perturbative regimes. The fac- 
torization formalism of nonrelativistic QCD (NRQCD) [1] provides a rigorous theoretical 
framework to describe the heavy-quarkonium production and decay by separating the am- 
plitude into two parts: a short-distance part which can be expanded as a sum of terms in the 
power of a s and calculated perturbatively, and long-distance matrix elements (LDMEs), 
which can be extracted from experiment. The relative importance of the LDMEs can be 
estimated by means of velocity scaling rules [2]. A crucial feature of the NRQCD is that 
the complete structure of the quarkonium Fock space has been explicitly considered. 

By introducing the color-octet mechanism (COM), NRQCD has successfully absorbed 
the infrared divergences into P-wave [1,3-6] and D-wave [7,8] decay widths of heavy quarko- 
nium, which can not be handled in the color-singlet mechanism (CSM). The COM can suc- 
cessfully reconcile the orders of magnitude discrepancies between the experimental data of 
J/tp production at the Tevatron [9] and the CSM theoretical predictions, even if it has been 
calculated up to the next-to- leading-order (NLO). Substantial progress has been achieved 
in the calculation of high order QCD corrections to charmonium production in order to 
clarify the validity and limitation of the NRQCD formulism. Recently, the complete NLO 
calculation for polarization and unpolarization of direct hadroproduction in NRQCD are 
presented by two groups [10-13]. For unpolarization of direct J /if) production in Ref. [11], 
two new linear combinations of color-octet matrix elements are obtained from the CDF 
data by the authors, and used to predict J/ip production at the LHC, which agree with the 
CMS data. J/t/j polarization puzzle also may be understood as the transverse components 
canceling between 3 S^ and 3 Pj 8 ^ channels [10]. The other group performed a multi-process 
fit of the color-octet LDMEs for unpolarization of direct hadroproduction [12], and gave 
more global fitting results which are consistent with NRQCD scaling rules. While, they 
find that the CDF polarization data (Run-II) [14] can not be interpreted by using these fit 
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values [13]. Therefore, the existence of the COM is still under doubt and far from being 
proven. 

The heavy quarkonium production associated with a gauge boson may be a good pro- 
cess in testing the COM. In high energy collider experiments, W, Z and heavy quarkonium 
can be selected with purely leptonic decay modes [15], which are particularly useful in 
hadron collides, because they provide an enormous suppression of the background. The W 
or Z production associated with a J/ip has been studied [16]. In reference [16] the authors 
gave a theoretical prediction for associated production of heavy quarkonim and electroweak 
bosons at hadron collider at the LO, and the numerical results show that the transverse 
momentum distribution of heavy quarkonium production is smaller in the CSM than that 
in the COM at the LO. Recently, the NLO QCD corrections to J/V> + W and J/ip + Z 
associated productions at the LHC were calculated in nonrelativistic QCD [17,18]. The 
numerical results show that the differential cross section at the LO is significantly enhanced 
by the NLO QCD corrections. 

Since the mass of bottomonium is about 3 times of that of charmonium, the value of 
v 2 is smaller in bottomonium (~ 0.1) than that in charmonium (~ 0.3), the expansion in 
v 2 should converge faster in bottomonium. Therefore, the study of bottomonium may be 
a more suitable choice to test the NRQCD factorization formalism. The Y(IS') production 
associated with a W- or a Z-boson at hadron colliders has been studied in Ref. [15]. It 
was found that the cross sections for T(IS') + W and T(15) + Z productions are roughly 
0.45 pb and 0.15 pb at the Tevatron and roughly 4 pb and 2 pb at the LHC respectively, and 
the dominant production mechanism is via the P-wave bottomonium production, which is 
a bb bounded state of color-octet, and the bottomonium subsequentially decays to T(IS'). 
In experimental aspect, this process has been studied once at the Tevatron [19], but due 
to the limited sensitivity of the CDF experiment in Run I, they didn't observe T(IS') 
mesons production in association with W- or Z-bosons at the Tevatron. Since the LHC 
has larger integrated luminosity and higher energy, we expect LHC experiments can achieve 
the sensitivity sufficient to observe the T(IS') plus a vector boson production. As we know, 
the NLO QCD corrections to quarkonium production are normally important in precision 
measurement [10-13]. It is necessary to provide a complete NLO theoretical prediction 
for T(IS') mesons in association with W- or Z-boson production, and compare the results 
with experiment. 

In this paper we calculate the full NLO QCD corrections to the T(IS') production in 
association with a VF-boson at the LHC. For this process, only the 3 Si color-octet (the 
bb 3 s[ 8 ^ Fock state) provides contribution at the leading-order (LO). Even including the 



NLO QCD contributions up to and v 7 of LDMEs, there are only color-octets bb 



lc-(S) 
^0 



do 



3 S^ and bb 3 Pj 8 ^ (J = 0, 1, 2), but no color-singlet contribute to the direct prompt 
T(15) production in association with a W gauge boson process. Therefore, the T(1S) + W 
production at the LHC can provide an ideal ground to study the COM. Meanwhile, The 
production of T(15) + W may also provide lampposts for looking for new physics [15]. In 
the context of supersymmetric models (SUSY), a charged Higgs boson lighter than about 
180 GeV have a sizable decay branching ratio into T(1S)W pairs if the ratio of vacuum 
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expectation values tan(/3) is small. Therefore, the T(IS') + W production may be a very 
important background process for searching for charged Higgs boson in SUSY models. 

The paper is organized as follows. We present the details of the calculation for the 
process pp — > T(1S) + W + X in Section 2. In Section 3, we give the numerical results and 
discussions at the LHC. Finally, a short summary is given. 



2. Calculation descriptions 

The strategy of the NLO QCD calculations in this work is similar to that used in the 
evaluations for the pp — > J/ip + W/Z process [17, 18]. For simplicity, we sketch it in this 
section. 

The cross section for the pp — > H + W process is expressed as 

, {PP ^H + W + X ) = J ^ £ N X) l Z(<>f " i[n] + W + X) 

x [G i/A (x 1 ,n f )G j/B (x 2 ,Hf) + (A^ B)] , (2.1) 



where N co i(n) and N po i(n) refer to the numbers of colors and polarization bb(n) states 
produced [6], respectively, a (i + j — > bb[n] + W + X) is the cross section for the short 
distance production of a bb quark pair in the color, spin and angular momentum state 
n at the parton level. (0 H [n]) is the long distance matrix element, which describes the 
hadronization of the bb quark pair into the observable bottomonium state H. G,ij/^^ are 
the parton distribution functions (PDFs), A and B refer to incoming protons at the LHC. 
i, j represents gluon and all possible light quarks and anti-quarks (i, j = u, d, s, c, u, d, s, c). 
We ignore the processes which are suppressed by a tiny Kobayashi-Maskawa factor. 

The short distance cross section for the production of a bb pair in a Fock state n, 
<r[i + j — > bb[n] + W + X], is calculated from the amplitudes which are obtained by applying 
certain projectors onto the usual QCD amplitudes for open bb production. In the notations 
of Ref. [6] they are written as 



Tr 



q=0 



A bhps{ 1/8) ] ~ £aTr 



c 1/8 n<?A 



9=0 



d 

'dq a 



C 1/8 U A 



9 =0 



A 



d 



Tr 



c (i/8) n i-4 



9=0 



(2.2) 



where A denotes the QCD amplitude with amputated bottom spinors, the lower index q 
represents the momentum of the bottom quark in the bb rest frame. n /i are spin projectors 
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Figure 1: Some representative Feynman diagrams for the process pp — > Y(IS) + W + X at the 
LO and up to QCD NLO. 



onto the spin singlet and spin triplet states. C\i% are color projectors onto the color singlet 
and color octet states. £ a and £ a p represent the polarization vector and tensor of the QQ 
states, respectively. 

In the calculation for the LO of a s , only the Fock state 3 s[ 8 ' 1 is involved. However, in 
considering the NLO QCD corrections to pp — > T(1S) + W + X process, we should not 
only consider the contribution from virtual corrections, but also the contribution from real 
gluon/light-quark emission processes. We draw some representative Feynman diagrams for 
the process pp T(1S) + W + X at the LO and up to QCD NLO in Fig.l. The virtual 
corrections only occur in connection with the bb Fock state 3 s[ 8 \ but the real gluon/light- 
quark emission process should involve 1 S , q 8 \ 3 S^ and 3 Pj 8 ^ Fock state contributions. In 
our calculations, we adopt the dimensional regularization (DR) method in D = 4 — 2e di- 
mensions to isolate the ultraviolet (UV), soft infrared (IR), and collinear IR singularities. 
There exist UV, soft, collinear and coulomb singularities in virtual corrections. In order to 
remove the UV divergences, we employ the modified minimal subtraction (MS) scheme to 
renormalize and eliminate UV divergency. After applying the renormalization procedure 
the UV divergences in the virtual correction part are canceled. The other appeared sin- 
gularities can also be analytically canceled when we added all the NLO QCD contribution 
components together. In Fig. 2, we present the IR and Coulomb singularity structures and 
divergence cancelation routes in the NLO QCD calculation for the pp — > Y(IS') + W + X 
process. 

The real gluon process with 1 S'q 8 - ) Fock state is free of divergence, while the contribution 



- 4 - 



Loop correction 


soft 


collinear 


coulomb 



Absorded in PDF 



finite 


soft colli near 


soft 




3 S (8) 


J 


Real gluon correction 




canceled by 
renormalized 
operator 

<o Y [ 3 s< 8) ]> 



-2L 



Absorded in PDF 



finite 


collinear 


finite 


1q(8) 


3g(8) 


3p(8) 


Real light-quark correction 



Figure 2: The IR and Coulomb singularity structures in the NLO QCD calculations for the 
pp -> 1(15) + W + X process. 



from the real gluon process with 3 S[ Fock state has soft and collinear singularities, and 
the real gluon processes with 3 Pj 8 ' > Fock states has soft singularities. The soft divergences 
from the one-loop diagrams will be canceled by similar singularities from the 3 s[ 8 ^ state 
contribution of soft real gluon emission. As for the collinear divergences from the 3 5j 8 ' ) 
state contribution of soft real gluon emission corrections, part of it can be eliminated 
by collinear singularities in virtual corrections, and the remaining collinear divergences 
in real gluon corrections can be absorbed into the PDFs. Nevertheless, it still contains 
coulomb singularities in virtual corrections and soft singularities arising from the 3 -Pj 8 ' ) 
state contribution of real gluon process, these singularities is neither collinear nor infrared 
in the usual sense, and it can only be eliminated in the spirit of the factorization approach 
by taking the corresponding corrections to the operator < O r ( ls ^[ 3 S^] > into account. 
As for the light-quark emission process, there only involve initial state collinear divergences 
in 3 s[ 8 ^ state subprocess, and it also can be absorbed into the PDFs. Finally, after adding 
all contributions together, a result which is UV, soft, collinear, Coulomb finite is obtained. 

We adopt the expressions in Ref. [20] to deal with the IR divergences in Feynman 
integral functions, and apply the expressions in Refs. [21-23] to implement the numerical 
evaluations for the IR safe parts of N-point integrals. In the virtual correction calculation, 
we find that only Fig. 1(e) and Fig. 1(f) induce coulomb singularities, and we use a small 
relative velocity v between b and b to regularize them [24]. The two cutoff phase space 
slicing method (TCPSS) [25] has been employed for dealing with the soft and collinear 
singularities in real gluon/light-quark emission corrections. 



3. Numerical results and discussion 



In the numerical calculations for the pp — > T(1S) + W + X process at the LHC, we 
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H 


< Of > < H [ 3 s} 8) ] > < H [ 1 5j 8J ] > 


T(1S) 
T(2S) 

Xo(lP) 
Xo(2P) 


15 x 1(T 2 GeV 3 2.0 x 1(T 2 GeV 3 
4.5 x 10~ 2 GeV 3 0.6 x 10~ 2 GeV 3 

2.03 GeV 3 4 x 10~ 2 GeV 3 

2.57 GeV 3 6.5 x 10" 2 GeV 3 



Table 1: NRQCD matrix elements for bottomonium production, where colour-octet matrix 
element < O h [ 1 S ( q S) ] >=< O H [ 3 P (8) ] > /m 2 b has been assumed for simplicity, and m b = 4.75 GeV. 

take CTEQ6L1 PDFs with an one-loop running a s in the LO calculation and CTEQ6M 
PDFs with a two- loop a s in the NLO calculation [26], and the corresponding fitted value 
a s (M|) = 0.130 and a s (M|) = 0.118 are used for LO and NLO calculations respectively. 
We take the number of active flavor as nj = 4. The relevant quark masses and fine structure 
constant are taken as: m q = 0, (q = u,d,c,s), mw = 80.398 GeV, = m T ^ ls ^/2 = 
4.75 GeV, a = 1/137.036. The renormalization, factorization, and NRQCD scales are 

chosen as /j, r = \if = mr and fi\ = mb, respectively, where tot(= yj (p^ 15 ^) 2 + 4m 2 ,) is the 
T(15) transverse mass. The NRQCD matrix elements of T(15) and other bottomoniums 
used in our calculations are list in Table [l], which determined from the CDF data [27,28]. In 
the calculation of real corrections, we adopt the two-cutoff phase space slicing method [25]. 
The two phase space cutoffs 5 S and 5 C are chosen as 6 S = 10 -3 and 5 C = 5 s /50 as default 
choice. In checking the independence of the final results on two cutoffs 5 S and 5 C , we find 
the invariance with the S s running from 10 -3 down to 10 -4 within the error control (less 
than one percent). Considering the validity of the NRQCD and perturbation method, we 
restrict our results to the domain p T > 3 GeV and |yx(is)l < 3. 

The dependence of the cross section on the renormalization scale \i T and factorization 
scale [if induces uncertainty for theoretical prediction. In Fig|| the \i dependence of the 
LO and the NLO QCD corrected cross sections with the constraints of p t > 3 GeV 
and |yx(i5)l < 3 for pp — > T(1S) + W + X process, is shown with our default choice fj, = 
fj, r = [Xf and the definition of fio = tut- There /i varies from //o/4 to 4/zq- Not like the 
usual expectation, Fig.^ shows that the NLO QCD correction can not improve the LO scale 
independence for the pp —> T(liS') + W + X process at the LHC. The related theoretical 
uncertainty varies from +8.1% to —8.6% at the LO and from +18.5% to —13.8% at the 
NLO when fi goes from ^o/2 to 2^o- Actually, this behavior is similar to that of the Wbb 
production at the LHC (see Ref. [29]). 

In our numerical calculations, we add all the contributions of pp — > T(1S) + W + + X 
and pp — > T(1S) + W~ + X together. In Figs.H(a) and (b), we present the LO and the 

T(1S) 

NLO QCD corrected distributions of p T and the corresponding K-factors for the direct 
T(15) + W production at the y/S = 8 TeV and the y/S = 14 TeV LHC, respectively. 
For comparison, we also depict the contributions from the 66[ 1 S'q 8 ' > ], 66[ 3 S'| 8 ^] and bb[ 3 Pj 8 ^] 
Fock states in these figures. From the figures we can see that the LO and the NLO 
corrected differential cross sections are sensitive to pt of T(IS'), and the LO differential 
cross section is significantly enhanced by the NLO QCD corrections, and the 66[ 3 5| 8 ^] 
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Figure 3: The dependence of the LO and the NLO QCD corrected cross sections for the process 
pp —> T(15) + W + X on the factorization scale and renormalization scale (/i/Mo) at the \/S = 
14 TeV LHC where we define fj, = jif = ji r and = m T- 




p; <iSl (GeV) p" ,s '(GeV) 



Figure 4: The LO and NLO QCD corrected distributions of p T and the corresponding K-factor 
for the direct T{IS) + W production, (a) at the \/S = 8 TeV LHC. (b) at the y/S = 14 TeV LHC. 



contribution dominates in the whole plotted p T region. In FigH(a), the NLO QCD 
corrected differential cross section of 1S * at the \fS = 8 TeV LHC decreases from 
0.126 pb/GeV to 0.0026 pb/GeV as p T increases from 3 to 50 GeV. In the range of 
3 GeV < p^ 1S ^ < 50 GeV, the K-factor, defined as K = rfg " Y (is) / d r^is) ; varies in the 

range of [2.85, 4.96], and reaches its maximum when p^ 1S ^ = 3 GeV. From Fig.|](b), 
we find that the NLO corrected differential cross section at the V~S = lATeV LHC in the 

Tf IS') 

whole plotted p T region can be quantitatively beyond 3 times of the corresponding LO 
differential cross section, and reaches its maximum K = 5.33 when p T = 3 GeV . We 
can read out from the figure that the NLO QCD corrected differential cross section varies 
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in the range of [0.211, 0.0053] pb/GeV whenp T runs from 3 GeV to 50 GeV, it reaches 

Tf IS) 

the maximum at the location of p T = 4 GeV. 

We have investigated the direct production of T(IjS') meson up to QCD NLO. The 
T(1S) meson can also be produced indirectly via radiative or hadronic decays of heavier 
bottomonium, such as T(2S), T(3S*), xm(1-P)> Xf>2(l-P)> Xf>i(2-P), and Xb2(2P) mesons. We 
know that these contributions are very important in studying T(IS') meson production [15]. 
To calculate these feed-down contributions for T(15) meson production, the main problem 
is to determine the relation between momentum of higher excited states and momentum 
of T(1S). Following the method used in Ref. [30], we ignored the momentum shift when 
these excited states decay into T(15). The cross sections of these six feed-down production 
channels can be obtained approximately by multiplying the direct-production cross section 
of the intermediate bottomonium production with its decay branching fraction to T(IS') 
meson. In our calculation, we ignored the contributions aroused from Xbo(^P,2P) feed- 
down into T(IS') for the smallness of the transition branching ratios. 

The direct-production cross sections of the T(2S) and T(3S') can be calculated by 
using analogous method to the T(15) production, and the relevant color-octet matrix 
elements of T(25) and T(35) used in our calculations are determined from the CDF data 
in Table |l] [27, 28]. In the calculation for the LO and the virtual corrections for the 
direct-production cross sections of the Xbii^-P)-, Xb2(^P), Xbii^P) and Xbi&P) mesons, 
only the Fock state is involved. However, in considering the real gluon/light-quark 

emission corrections process, we should not only consider the contribution from 3 5ij Fock 
state contribution, but also the contribution from 3 P^ Fock state contributions. Based 
on NRQCD factorization formalism [6,31], we can handle these calculations by adopting 
analogous method to the direct T(IS') production. For the NRQCD matrix elements 
of Xbj, the color-singlet matrix elements are taken from the potential model calculation 
of Refs. [27,32,33], and the color-octet matrix elements are determined from the CDF 
data [27,28]. The branching ratios of heavy bottomonium decays to T(1S) mesons used in 
this paper are taken form Ref. [34]. 




Y (IS) 

Figure 5: The complete NLO QCD corrected distributions of p T for the pp — > T(IS') + W + 
X process, (a) at the y/S = 8 TeV LHC. (b) at the y/S = 14 TeV LHC. 

The complete NLO QCD corrected distribution of for the pp — > T(IS') + 

W + X process, the contributions of the T(IS') direct-production and the T(1S) indirect- 
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production by T(2S), T(3S), x&(l-P), and Xb(2P) mesons feed-down decays at the LHC 
are illustrated in Fig.[|. We can see that the T(15) indirect-production can give important 
contribution to the distribution of p T for the pp — > T(1S) + W + X process at the 
NLO. We find that when px is smaller than about 7 GeV , the contribution from Xb(l-P) 
feed-down decays is about the same order of magnitude with the contributions of T(1S) 
direct-production when \/S = 8 TeV, particularly at the V^S* = 14 TeV LHC the produc- 
tions of Xb(lP) feed-down decays provide the main contribution to p T distribution for 

T( IS) 

the pp — > T(1S) + W + X process. With the increment of the p T , the contributions of 
T(1S) indirect-production contributions decrease quickly, and the T(1S) direct-production 
contribution dominates at the large pt region. 

4. Summary 

In this paper we investigate the complete NLO QCD correction of T(15) + W associated 
production at the LHC. This process is an ideal platform for studying color-octet mecha- 
nism. We adopt the dimensional regularization to deal with the UV and IR singularities in 
our calculation. The Coulomb and soft singularities in P state are isolated and absorbed 
into the NRQCD NLO corrected operator < O r( - ls ^[ 3 s[ 8 '] >. After adding all contribu- 
tion components together, we get the results with UV, IR, Coulomb safety. We find that 
the production rate of the pp — > T(15) + W + X process is quite large, and this process 
has the potential to be detected at the LHC. For the T(15) direct-production, the dif- 
ferential cross sections at the LO are significantly enhanced by the QCD corrections, and 
the bb [ 3 S\ ] contribution dominates in the range of 3 GeV < pp 1 < 50 GeV. In this 
paper, we also calculate the T(IS') meson indirect-production via radiative or hadronic 
decays of T(2S), T(35), X6l(LP), Xm(1P), X&i(2P), and xm(2P) mesons. We find that 
the T(IS') indirect-production contributions can give important contribution to the distri- 

T( IS) 

bution of p T for the pp — > T(IS') + W + X process at the NLO. In the lower px region, 

T(1S) 

the contribution from indirect-production can give the main contribution to p T in the 
pp — > T(15) + W + X process. We conclude that if the T(IS') + W production is really 
detected at the LHC, it will be useful for testing the NRQCD factorization formalism. 
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